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SUMMARY 


Resistance  thermometer  corralation  measurements  o £  the  thermal 
wakes  behind  two  line  heat  sources  are  suggested  for  estimation  of  the 
(Lagranglan)  Joint  probability  density  of  fluid  particle  displacements.  Con¬ 
siderable  uncertainty  results  in  practice  from  the  concomitant  molecular 
heat  diffusion. 

Experiments  Indicate  that  in  addition  to  the  well-known  nearly  normal 
(Oaossian)  shape  of  the  probability  densities  of  individual  particle  displace¬ 
ments*  these  displacements  are  also  close  to  being  jointly  normal. 
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INTRODUCTION 

If  we  mark  with  dye  a  particular  portion  of  a  turbulent  flow,  we  shall 
••a  the  marked  fluid  change  position  and  shape  as  time  progresses.  These 
two  changes  are  manifestations  of  the  diffusive  character  of  the  flow.  If  we 
repeat  the  ejqpariment  in  a  "statistically  Identical"  flow,  the  changes  expert* 
enced  by  the  marked  fluid  will  in  general  be  different,  but  if  the  experiment  is 
repeated  often  enough  certain  average  quantities  can  be  defined.  The  same 
will  be  true  in  a  single  flow  with  statistical  properties  stationary  in  time.  For 
example,  the  mean  square  displacement  of  the  fluid  mass  from  the  mean  path 
line  and  the  average^ increase  of  a  dimension  of  the  marked  fluid  volume  are 
,  two  such  quantities  characterising  the  diffusion.  It  is  the  goal  of  turbulent 
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diffusion  theory  and  experiment  to  be  able  to  predict  such  (Lagrangian)  measures 
*  of  diffusion  from  other  properties  of  the  flow  field,  preferably  those  more  easily 

i 

measured. 

When  the  contaminant  doen  not  change  the  dynamical  properties  of  the 
marked  fluid,  this  problem  is  obviously  kinematlcal  (in  the  absence  of  molecular 
transport):  knowledge  of  the  properties  of  path  lines  in  the  fluid  would  be  sufficient 
to  determine  the  turbulent  convection. 

Taylor^  suggested  that  even  with  molecular  transport,  the  obvious  in- 
'  dependence  of  random  molecular  migration  and  random  turbulent  migration 
permits  simple  superposition  of  the  two  resulting  mean  square  displacements, 
hence  the  two  transport  effects.  Townsend'2)  has  recently  disproved  this  con¬ 
clusion  on  the  grounds  that  the  obvious  distorting  effect  of  turbulent  convection 
on  local  scalar  gradients  actually  changes  the  local  molecular  transport  rate. 


The  problem  of  diffusion  in  a  homogeneous  turbulent  flow  with  no 
molecular  diffusion  was  treated  by  Taylor  in  1921^,  He  related  the  simplest 
measure,  the  mean  square  displacement  of  a  single  fluid  particle  as  a  function 
of  time,  to  the  double  Lagrangian  time  correlation  of  velocity.  So  far  no  one 
has  been  successful  in  predicting  this  correlation  function  from  the  general 
dynamics  of  the  flow  or  in  relating  it  to  the  (more  accessible)  Eulerian  properties 
Its  traits  are  known  only  through  experiment.  Another  pertinent  function,  the 
probability  density  of  particle  displacement  in  a  homogeneous  turbulence,  has 
been  shown  experimentally  to  be  normal  at  all  diffusion  times^'  but  no 
theory  has  been  able  to  deduce  this  fact. 

At  the  present  time,  therefore,  we  have  available  Taylor's  theoretical 
connection  between  diffusion  and  Lagrangian  correlation  [written  in  more  general 
form  by  Batchelor^  J  plus  the  empirical  evidence  for  a  normal  probability 
density  of  particle  displacement  in  shear-free  homogeneous  turbulence.  By 
superposition  of  point  source  results  these  are  sufficient  for  prediction  of  the 
average  contaminant  concentration  field  due  to  any  prescribed  spatio-temporal 
distribution  of  source  strength. 

Taylor's  approach  is  applicable  only  in  a  homogeneous  turbulence,  and 

/ 

the  most  detailed  data  are  available  in  the  approximately  homogeneous  (and 
isotropic)  turbulence  behind  a  regular  grid  in  a  wind  tunnel.  Our  resulting  semi- 
empttical  understanding  supplies  us  at  least  with  a  jumping  off  point  for  conjecture 

and  for  inspection  of  the  diffusion  In  more  complex  "practical"  turbulent  flows*, 

■  ■  - ■■  —  -  ■  "■  "  - - — .  . .  . 

*See,  for  example,  "Remarks  on  Turbulent  Heat  Transfer"  by  S.  Corrsin, 

Proc.  Iowa  Conference  on  Thermodynamics,  April  1953. 


When  more  detailed  information  is  desired,  for  example  the  probability 
that  the  concentration  at  a  fixed  space  point  will  exceed  a  chosen  value  with 
more  then  one  active  source,  superposed  single  source  data  are  inadequate;  we 
must  know  the  Joint  probability  function  for  the  particle  displacements  from  the 

several  sources. 

The  simplest  generalisation  of  the  fairly  well -studied  single  line  source 
case  is  that  of  two  parallel  line  sources.  The  simplest  turbulence  available  is 
the  decaying  approximately  isotropic  turbulence  behind  a  grid.  The  most  con¬ 
venient  contaminant  is  heat.  The  present  work  was  undertaken  to  establish 

an  experimental  method  for  measuring  the  joint  probability  function  for  such 

\ 

an  arrangement. 

This  Joint  dispersion  problem  has  been  analytically  described  in  some 
detail  by  Batchelor^81  who  also  treated  special  aspects  theoretically^),  In  both 
works  he  excluded  molecular  diffusion  effects.  The  paper  of  Brier^)  should 
also  be  noted. 

The  problem  and  experimental  technique  developed  here  were  suggested 
by  Professor  S.  Corrsin,  who  also  supervised  the  investigation.  Miss  V.  O'Brien 
carried  out  much  of  the  calculation. 


THE -PKQBLEM  OF  THE  LINE  SOURCE  IN  TURF^ENT  FLOW  WITH 

MOLECULAR  DIFFUSION 

With  the  usual  average  space -time  transformation  used  in  isotropic 

•  the  spread  of  the  region 
contaminated  by  the  marked  particles  with  increasing  downstream  distance 
x'rom  the  source  is  a  measure  of  their  lateral  diffusion.  This  diffusion  is  not 
exclusively  a  result  of  the  turbulence,  however,  since  even  if  the  turbulence 
level  were  sero  there  would  be  a  spread  due  to  molecular  effects.  If  there  were 
no  moleciilar-iliffuBion,  the  marked  particles  would  constitute  a  surface  with 
wrinkles  which  tend  to  increase  with  time,  hence  with  downstream  distance.  The 
general  spatial  region  occupied  intermittently  by  this  wavey  sheet  we  call  the 
"turbulent  wake".  The  local  region  surrounding  the  particles  marked  at  the 
source,  where  the  "marking"  diffuses  by  molecular  motions^  will  be  called  the 
"molecular  wake";  hence  this  is  the  Instantaneous  wrinkled  sheet  as  thickened 
by  molecular  transport.  The  total  region  influenced  by  the  "marking",  irre¬ 
spective  of  the  diffusive  mechanism,  will  be  called  the  "mean  wake". 

The  word  "wake"  here  lo  used  to  denote  the  region  contaminated  by  the 
source.  The  fluid  particles  defining  this  region  are  marked  in  come  detectable 
fashion,  in  this  case  with  heat.  The  sources  and  the  contaminant  are  ansumedto 
have  negligible  effect  on  the  dynamical  properties  of  the  fluid,  as  restriction 
satisfied  by  the  actual  experiments;  the  momentum  wake  was  undetectable  \f?.  inch 
behind  the  source  wire. 

In  accordance  with  a  suggestion  of  Taylor^,  the  two  mechanisms  causing 
diffusion  are  usually  asaumed  to  be  uccorrclated  in  the  sense  that  the  mean  square 
width  of  the  mean  wake  is  taken  to  be  equal  to  the  cum  of  the  mean  square  widths 

of  the  turbulent  and  the  undisturbed  molecular  wakes.  Townsend^  has  recently 
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pointed  out  that  they  are  actually  correlated  in  the  following  sense:  the  molecular 


turbulence  experiments 


i-  ■==.  u 

dLt:  ckx- 


di/furion  in  not  independent  of  the.  urbiv  on1;  rno  iin  because  .'u  .w-sn  turbulent 

•train  convecllvely  change  the  instantaneous  temperature  gradients  in  such  &  way 
that  the  average  molecular  wake  width  a£;er  a  given  diffusion  time  is  greater  than 
that  of  the  molecular  wake  in  a  non-turbulcnt  flow. 

If,  however,  it  is  found  that  the  actual  average  molecular  wake  width  is 
the  .same  regardless  of  its  centerline  location  (i,  e. .  of  its  position  within  the 
turbulent  wake),  then  Taylor's  assumption  can  he  reinvoked  in  modified  form: 
ths  turbulent  dispersion  is  uncorrolatcd  with  the  accelerated  molecular  dispersion. 
Henee  the  squares  of  the  corresponding  otandar  i  deviations  i  add  to  give  the 
equated  standard  deviation  of  the  mean  wake: 


In  Taylor's  original  suggestion^)  0^  vas  taken  from  the  non-turbulent 
molecular  wake. 

In  this  investigation  the  above  lack  of  correlation  is  established  experi¬ 
mentally.  A  theoretical  proof  would  involve  shoving  zero  correlation  between 

r* 

the  turbulent  particle  displacement  J  ~  j  V*-  (  i , )  cit , 

«> 

and  whatever  function  of" strain  rate  r \j  velocity  spatial  derivative  is 

proportional  to  the  increase  in  average  molccuJar  wake  width. 

Townsend's  accelerated  molecular  diCfu  dor  as  described  above  is  easily 
estimated  only  for  turbulence  spectral  componc  its  (eddies  V)  with  wave  length 
much  larger  than  the  molecular  wake  width,  T'  is  per-nit3  p.ppreximatior.  by  r. 
spatially,  constant  strjxiELfleld*.  For  spectral  components  with  wave  lorphh 


*  See  the  analysis  in  Appendix  2  of  reference  Z, 


much  smaller  than  the  molecular  wake  a  different  viewpoint  scemn  ?pprppriate: 
these  fluctuations  spread  the  molecular  wake  from  "within"  just  as  the  molecular 
motion  does.  Then  the  term  "molecular  wake"  loses  whatever  appropriateness 
may  have  remained.  The  effective  diffusivity  of  this  very  small  scale  motion 
could  be  estimated  with  the  analysis  of  Corrain^*)  for  heat  transfer  in  isotropic 
turbulence  with  prescribed  mean  temperature  gradient. 

For  turbulent  fluctuations  of  scale  comparable  with  the  molecular  wake 
width  neither  of  the  above  simplified  accounts  is  adequate,  and  a  full-fledged 
non-uniform  strain  analysis  would  have  to  be  carried  out. 


ANALYTICAL  .REPRESENTATION 


'  1 


Taylor's  now  classical  analysis  of  the  diffusion  from  \  single  line  source 
in  homogeneous'  non-decaying  turbulence  with  no  molecular  diffusion  shows  that 
the  mean  square  spread  of  the  turbulent  wake  is  given  by^) 

i  II!  s  2  7*  \S(  t  )  *  zF  fJESDplIi^  ar 

?  at  J  i  v*- 

?  ___  * 

k 

where  \j  *  fr(k,)Ji is  particle  displacement.  \T  la  fluctuation  velocity 

t  o 

in  the  y  -direction. 

|  Perhaps  the  main  problem  in  single  source  diffusion  theory  is  the  prediction 

j*!.  ' 

of  f  £  (imt)  •  tho  Lagrnngian  time  correlation,  in  terms  of  the  more 

$ 

i  S 

easily  measured  Eulertan  statistical  function*.  The  analysts  leading  to  the  above 

I 

requires  the  turbulence  to  be  stationary  in  time.  Therefore  diffusion  in  decaying 

I 

turbulence  can  be  fitted  into  thlt,  simple  mold  only  if  it  obeys  some 

suitable  similarity  conditions.  Townsend^  has,  in  effect,  assumed  complete 

| 

similarity  of  all  functions  both  Eulerian  and  Lagrnngian.  This  assumption 

I 

collapses  the  diffusion  data  remarkably  well  conr.1  dering  that  there  r-eemo  to  be 

\ 

no  a  priori  reason  to  expect  such  detailed  self-proBcrvation  of  the  do  raying 

turbulence  field. 

|  This  gives  a  result  for  $  ft)  approximately  independent  of  the  locution 

f 

of  the  source  In  a  given  decaying  field.  Before  the  calculation  can  bo  carried  cut 

t  — 

the  molecular  Influence  has  to  be  removed  from  the  y  A  data,  T’hich  he  did 
by  measuring  the  accelerated  tnoleculrr  '  iifunton  and  by  iioaumtng  i  ;  uncorratntad 

k 

i  - 

to  the  turbulent  diffusion. 


{ 

If  the  functional  form  of  the  probability  density  of  the  displace; -tent  ic 

characterised^byjane  parameter  (such  as  its  standard  deviation)  then  a  knowledge 

1  . 

of  the  Lagranglan  correlation  function  gives  a  fairly  complete  statistical  picture 
of  the  diffusion.  In  "isotropic"  wind  tunnel  turbulence  this  probability  density 
Auction  is  observed  to  be  normal  (Gaussian)  at  all  diffusion  times  so  that  know- 

i 

ledge  of  the  Lagrangian  correlation  function  permits  the  mean  concentration  to 

•  . ,  i  '  a 

be  computed  at  a  given  point  for  any  source  distribution. 

It  is  pertinent  to  examine  the  timo  variation  of  temperature  at  a  fixed 

«■ 

point  in  the  wake.  Near  the  source  this  time  variation  is  very  spiky  with  a  fixed 
variation  between  the  ambient  temperature  (say  zero)  and  some  peak  value.  ^ 

X 

This  peak  is  fairly  insensitive  to  velocity  gradients,  so  that  any  observed  variation 
4  < 

ctn  peak  temperature  at  the  fixed  point  is  caused  by  the  Molecular  wakes  not  being 

4 

Completely  swept  past  the  point.  Near  the  edge  of  the  turbulent  wake  the  tempera¬ 
ture  is  aero  except  for  infre'ment  pulses  due  to  the  molecular  wake's  entering 

this  region.  Near  the  center  of  the  turbulent  wake  the  pulses  are  more  frequent 

1 

and  the  temperature  only  rarely  drops  to  ths  zero  value.  After  greater  diffusion 
time  the  signal  has  essentially  the  same  character,  but  now  there  io  no  clearly 

defined  maximum  temperature  since  the  instantaneous  peak  value  is  determined 

.< 

by  the  tlmo  history  of  the  turbulent  and  laminar  effects  on  the  dlfl‘  usioa,  and 

/  \ 

the  longer  the  time  the  larger  the  variability.  For  example,  after  great  diffusion 

time  the  skewness  of  ths  temperature  fluctuation  at  a  fixed  space  point  is  con- 

|  - 

Aldsrably  reduced^). ^ 

4  The  source  configuration  considered  in  this  paper  Is  that  of  two  parallel 

i 

line  sources  at  the  same  downstream  distance  from  the  turbulence  producing 
grid.  The  problem  has  two  parameters  in  the  sense  that  the  spacing  of  the  sources 
and  the  distance  downstream  from  the  sources  are  the  variables  used  in  exploring 
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the  Joint  statistical  properties  of  the  wakes.  For  both  aero  source  facing  and 
•  large  source  spacing  the  problem  reduces  to  an  application  of  single  source 
results,  since  for  apactngs  large  compared  to  all  length  scales  of  the  turbulence 
the  wakes  move  independently. 

The  function  of  interest  In  this  investigation  is  the  joint  probability  density 
of  particle  displacements  from  the  two  sources.  As  will  be  eihown  shortly,  this 
is  approximately  proportional  to  a  measurable  quantity,  the  correlation  between 
temperature  fluction  signals  at  two  points  in  the  turbulent  wakes.  This  relationship 
is  the  essence  of  the  experimental  technique  introduced  here. 


In  particular r  if  the  stream  direction  is  X  ,  the  direction  perpendicular 
to  the  main  stream  and  the  line  sources  is  ,  and  the  direction  parallel  to  the 

sources  is  a  ,  then  the  function  of  interest  is  the  probability,  defined  as  a 
time  average  that  one  turbulent  wake  is  contained  in  a  region  (*,  y,,*)** *) 
and  simultaneously  the  other  turuulent  wake  is  contained  In  a  region 


(*  (*,Jt *  )  (rl«*  Z) '  thC  experlm0nt*  th*  *  -location 

of  the  sources  is  constant,  say  zero,  and  all  measurements  are  made  at.  the  same 

%  ,  Therefore,  thejprobability  Is  a  function  of  the  source  spacing  Y"  , 

the  distance  from  the  source  X  ,  and  the  two  lateral  positions, 


This  function  is  denoted  by 


P^.JzSY.x) 


where  \j  and  are  the  individual  ^  -displacements  of  the 

two  wakes. 

3y  a  priori  reasoning  it  is  possible  to  make  certain  statements  about  this 
function.  Near  the  source  the  probability  density  function  of  the  displacements 


) 


Is  Identical  to  that  of  the  velocities  at  :ie  two  pokta  :  we  r>p,\.y  the  :  p  •  >;prUt<; 

space-time  transformation.  Therefore,  all  known  statistical  properties  should 
be  recoverable  from  the  displacement  measurements,  e.g,  the  double  correlation 
and  the  skewness  of  the  velocity  difference.  Within  a  reasonable  distance  from 
the  sourcts,  for  a  fixed  3r  ,  the  two  wakes  cannot  cross,  and  hence, 

>  V,*)  a  0  kr  yx-j,  >  V 

The  traneveree  isotropy  of  the  flow  also  imposes  the  condition  that 

y,*) --  ;  y,  *) 

The  remaining  restriction  is  that  the  single  wake  should  be  recoverable  from 
the  marginal  distribution.  if  the  single  wake  probability  density  for  the 
displacement  is  |j  (y  j  %)  ,  than 

■i 

\ 

V 

f 

and 


b (»>)*)  *  I  P 

—  oo 

P  k»;x)  -  /  P  d.Y 

-  »o 
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EXPERIMENTAL  EQUIPMENT 


The  experimental  work  war  carried  out  in  a  2  x  2  foot  N.  P.L,  type 
open  return  wind  tunnel  at  a  mean  velocity  of  15  feet  per  second  (Pig.  1).  The 
turbulence  waa  generated  by  a  bl  -plane  square  me  ah  grid  with  1-inch  me  ah  and 
1/4 -inch  round  rods* 

Two  spring  loaded*  vertical*  .003 -inch*  heated  Nlchrome  wires  which 

spanned  the  tunnel  served  as  sources  for  marking  the  fluid  particles.  The  wires 

were  mounted  20  mesh  lengths  from  the  grid  in  a  frame  which  permitted  independent 

parallel  motion  in  the  vj  -direction.  The  wires  were  heated  with  drect  current 

^  \ 

to  a  temperature  of  approximately  300°  C*  this  temperature  producing  no  measurable 
dynamical  effects  on  the  heated  fluid  at  the  position  where  measurements  were 
taken*  There  was  no  observable  sag  or  motion  of  the  wires  due  to  aerodynamic 

forces* 

? 

[  The  instantaneous  temperature  was  detected  downstream  of  the  sources  by 
two  platinum  resistance  thermometers  ,00005  inch  in  diameter*  constructed  from 
Wollaston  wire.  The  thermometers  were  operated  with  sufficiently  low  current 
(  <  .,001  amp)  ao  that  the  velocity  output  was  negligible  compared  to  the  tempera¬ 
ture  output  at  the  X  -locations  used^^.  These  thermometers  were  mounted  on 
a  traversing  mechanism  that  permitted  independent  motion  in  the  -direction* 
i,  e.  perpendicular  to  the  mean  flow  direction  and  to  the  sources.  The  -location 

with  respect  to  the  source  location  could  be  established  within  .  002  inch. 

% 

♦ 

„  Mean  temperatures  were  measured  by  averaging  the  temperature  signal  with 
long  period  galvanometer.  One  thermometer  was  used  in  the  wake  system*  and 


11’. 

the  other  was  placed  in  the  stream  outside  the  region  oJ.'  wake  influence.  Large 
scale  temperature  fluctuations  affecting  both  wires  were  cancelled  by  a  bridge 
circuit.  A  Leeds  and  Northrup  K-2  potentiometer  was  used  to  measure  the 
voltage. 

The  instantaneous  wake  temperatures  were  converted  to  usable  voltages 
by  amplifying  and  compensating  the  thermometer  output  in  accordance  with 
ordinary  hot-wire  anemometry  practice.  The  thermometer  had  a  time  constant 
of  •  I  millisecond  audjthencompensated  voltage  output  was  flat  to  within  2%  from 
2  to  8000  cycles  per  second. 


-MEASUREMENT  OS  THE  JOINT  PROBABILITY  OF  '.'Hr.  V.'Afe. 

DISPLACEMENTS 

7 

The  joint  probability  of  partiele  (Replacement!  from  the  two  sources  i  t 
essentially  the  fraction  of  total  time  that  the  wakee  are  rfmultaneouily  at  two 
•pacified  v  •locations.  Since  our  thermal  tagging  diffuses  molecularly, 
correction  must  be.  made  if  we  seek  the  turbulent  dispersion  only.  Toward 

j>* 

this  end  we  need  a  measurement  of  the  average  "thickness”  of  this  hot  sheet 
for  eaeh  position  in  the  X  ,  y  plane.  In  fact  we  should  like  an  avarnpe 
''local'*  temperature  distribution  for  each  such  position.  This  can  bo  measured 
approximately  in  a  olnglo  wake  via  the  correlation  between  temperat-.iruafluct  latipnei 
picked  up  by  two  resistance  thermometers  (cool  "hot-wireo")  separated  hi  thr 

if 

y  •direction. 

I 

%  Let  £(^1  j  bo  the  single -wake  temperature  Held.  The  sen  i- 
colon  ia  simply  a  reminder  that  we  foous  attention  upon  r  (j)  for  fined  t 


and  ,  -ffy)  can  be  described  also  as  a  function  of  tlia  i  -location 

i  tne«.n 

of  Its  mean  (  )  and  the  distance  from  this^  (  y  -  JJf  )  ,  i.  e. 

•V-  —  ■  * 

I  ftj).  *  f,  (sr. r's  *,*) 


where 


./  s 


J  y  PtyMy 


* 

■  f  / 

|  For  each  value  of  y  (anc  x  of  courts)  tfe^ro  is  cn  r.'-emge. 

t  ' 


temperature  cietribution _ jp;  *  )  l  this 


ill  the  one  wo  ;jhcu  d 


like  to  determine  experimentally  for  the  purpose  of  correcting  the  e  :oca  mes  * 

wake.  In  principle  it  could  be  determined  as  a  time  average  of 


$1  ( *  )  y  *t  yf  t  J  with  j/  heldflxed,  l,e.  the  average  of 

the  collection  of  instantaneous  f 's  which  occur  whenever  the  mean  position 
of  the  hot  sheet  sweeps  past  a  chosen  y  f 

- - -  /V 

In  order  to  make  an  experimental  estimate  of  r;  we  assume 
(a)  that  -C,  is  a  normal  (Gaussian)  curve,  though  lower  and  broader  than  the 
normal  thermal  wake  which  would  occur  with  no  turbulence;  (b)  that  the  instanta¬ 
neous  p  (y)  profiles  differ  only  a  small  amount  from  the  "conditional 
mean  function"  *f t  ^  "  y  /  .  Under  these  assumptions,  the  directly 

measurable  temperature  correlation  function ^ 

is  approximately  equal  to  the  hypothetical  correlation  that  would  result  from 

**  . 

random  lateral  translation  (by  variation  of  y  *  )  of  a  fixed  profile  f(  (y-y'J 

OO 

«  f  r-y'J 

^  oO 

where  /*  (yO  is  the  probability  density  of  y  /  i„  e.  of  the  turbulent 


wake. 


With  assumption  (a)  plus  the  known  closely  normal  character  o>.  ( y  ') 


we  get 


*fr*‘*a^  ]&, 


/*  -  3  'J 1  (y  «;  <*_  -y ■> 1  y ' 

7.  er  *  “  TcrT"”  a  cr  »■  . 

e  '  .  e.  '  '  e  * 


OMrl  wf 


a  »  t  «  « 


r)  ^  Sv~«»)J  f 


where 


and 


^  are  the  standard  deviations  o£  f  auu  p 

What  we  want  to  check  te  whether  is  independent  of  J 

,<2) 


are  the  standard  deviations  of 


auu 


as  assumed  by  Townsend'  Since  the  independence  assumption  works  fairly 
well  for  single  wake  measurements, 

Anction  of  \j 

distance  comparable  with  cr  ,  whence  th~  major  contribution  to  the  integral 
comes,  then  we  can  write  approximately 


0*^  can  be  at  most  a  slowly  varying 
.  In  particular,  if  the  variation  of  cr  is  small  over  a 

% 


proximately  - 

K-5?  £> * '*  -  r"  ~  } 

R<',r)5S  ^<r  u 


Presuming  that  the  measured  •  ^  ^  r  i  .Y  )  does  give  normal  R^r) 

for  each  ^  ,  and  that  is  known,  th-  above  expression  could  be  need 

in  principle  to  calculate  ^  )  .In  fact,  Cy  is  not  known  a  priori 

but  we  do  have  <T  ,  the  standard  deviation  of  the  mean  thermal  wake.  For 
© 

this  calculation  we  must  therefore  estimate  0^  ;  the  simplest  estimate  is 

perhaps 


\ 


where  ^Crl\  is  the  average  value  of  (v)  acroa;.  the 

r  'Have ,  r  J 

mean  wake,  at  fixed —  X 

A  straightforward  way  to  use  this  statistical  information  on  -p  for 
determining  *P  Clfij  y»,  j  V*  %  )  is  to  measure  the  double  correlation  be¬ 
tween  the  temperature  fluctuations  in  the  two  walces,  the  measurements  being 
taken  at  the  same  *X  ,  The  result  of  such  a  measurement  would  be  a 
function  of  ,  the  thermometer  positions.  This  function  is 


5 


a',*)  'Jo v,*)  Jx/r,a*' 


This  expression  assumes  that  la  Independent  of  at  the 


aama  inatant 


.  ?  t. 


defined  aa  an  average,  and  thla  assumption  means 


that  the  average  ahape  of  -T(  at  a  particula  r-  y  la  independent  of  a 
aubaidlary  condition  auch  aa  the  fact  that  nearly  fluid  particles  have  also  been 


tagged. 


Unfortunately,  inveraion  of  the  above  expression  for 


■P  is  a 


imple  only 


for  special  easea,  auch  as  P,  being  a  Dirac  type  function;  hence,  thie 
approach  le  not  advisable  if  others  are  available.  Furthermore,  the  standard 
deviations  of  f*,  and  Jo  in  the  experimental  configuration  used  happened 
to  be  of  the  same  order  of  magnitude,  which  makes  far  from  a  Dirac 

function. 

Since  the  information  sought  requires  measurement  of  the  location  of  the 
fluid  which  has  actually  been  "in  contact"  with  the  sources,  a  method  was  devised 
that  used  only  the  peak  values  of  the  temperature.  The  turbulence  keeps  the 
molecular  wake  symmetric  on  the  average,  so  that  the  peak  temperature  is  located 
at  the  symmetry  point,  i.e.  the  point  that  was  "la  contact"  with  the  sor.ee. 

The  manner  in  which  the  measurements  were  performed  van  to  make  the 

* 

temperature  at  a  point  actuate  a  trigger  circuit  -  \  such  a  way  that  when  the  tempera* 
tore  exceeded  a  certain  level,,  a  constant  non-zr-o  voltage  was  produced  which  went 
again  to  zero  aa  soon  aa  the  temperature  clroppe  i  below  this  level,  V/lth  two 
channels  of  this  type,  one  for  each  wake,  signal*  are  generated  which  yield  a  fixed 
voltage  when  either  transducer  is  in  a  region  where  the  wake  exceeds;  a  given 
temperature.  These  binary  signals,  fluctuating  jetween  zero  and  a  fixed  value 


not  sero,  depending  on  the  absence  or  presence  of  the  hot  region?  c '  t!!.c  wakes, 
when  multiplied  together  yield  a  signal  whose  integral  in  proportional  to  the 
desired  Joint  density. 

This  type  of  operation  on  the  temperature  signal  converts  each  molecular 
wake  into  a  rectangular  shape  whose  width  £  is  determined  by  the  width  of 
the  molecular  wake  at  the  temperature  used  to  actuate  the  trigger.  If  this  rect¬ 
angle  Is  called  Q  \  if  j  •  the  number  obtained  by  multiplying  to¬ 

gether  tie  signals  from  both  wakes  is  given  by 


— - •• 

Si  i \*)  ~  JjV •Ss.'^sL 


If  the  critical  temperature  is  high  enough  so  thAt  the  variation  of  P  is  small 
in  the  region  £  ,  than 

Gj  -  v, ’•)«,«.  . . . 

Using  the  signals^in  the  manner  indicated  here  requires  that  an  accurate 
reproduction  of  the  temperature  signals  be  available  to  actuate  he  trigger.  The 
hot-wire  amplifiers  were  a.c.  coupled  for  reasons  of  stability,  so  that  the  out¬ 
put  signal  from  tho  amplifier  had  a  aero  mean  value.  In  order  to  regain  a  vol’&ge 
actually  proportional  to  the  instantaneous  temperature,  it  is  necessary  to  regain  the 
moan  value.  This  was  done  by  either  clamping  the  signal  t.o  the  ambient  tempera¬ 
ture  level  of  the  signal,  i.  e.  making  sero  voltage  correspond  to  the  ambient  tunnel 
temperature,  in  locations  where  a  well-defined  r.ero  er.ieted  (near  the  edges),  or 
by  adding  a  voltage  proportional  to  the  mean  temperature  in  regions  where  no 
bate  level  existed  (wske  center). 


r 


These  signal*  were  shaped,  as  described  above,  v/lth  Schmitt  V.:i;jger 
circuits  into  on -off  signals.  The  trigger  signaln  from  each  wake  were  instant¬ 
aneously  multiplied  and  the  resulting  coincidence  signal  was  integrated  to  give 
a  measure  of  the  joint  probability  distribution.  The  integration  was  accomplished 
by  using  the  coincidence  signal  to  modulate  a  regular  pulse  train  and  counting 
those  pulses  that  passed  when  the  signal  was  non-zero.  A  diagram  of  the  circuit 
that  performs  the  necessary  operation  is  shown  in  figure  (3). 

In  order  to  determine  the  correct  voltage  to  be  added  to  the  signal  for 

r 

restoring  the  mean  temperature  value,  measurements  of  the  wake  of  a  single 
line  source  were  made  with  the  same  method  used  for  the  two  wakes,  except 

t* 

that  hera  the  integration  was  performed  on  instead  of  the  coincidence 

signal  from  the  two  wakea.  The  result  of  this  measurement  Is  equivalent  to 


H's)  *  f  °cs’) 


M s'> 


was  obtained  from  the  mean  wake  and  the  orperimentally  determined 


4 l  •  Sufficient  voltage  wao  added  at  each  point  of  the  signal  to  make 

H^j)  proportional  to  •  In  general,  appreciable  acldll’on  of  voltage 

waa  needed  in  only  the  central  portion  of  the  turbulent  wake. 


RESULTS 


»  > 


Munr«n«8ti  were  mad*  at  tha  distance*  of  X  ■  7/8  ,  6-7/8,  and 
18*7/8  Inch**  from  tha  aourcaa. 

C  Tha  function  ffy  +  r)  was  measured  in  order  to  obtain  information 

about  tha  molecular  wake*  necessary  to  adjust  the  equipment  and  to  interpret  the 
results.  The  measured  points,  normalised  on  4s  tel  ?<e)  ,  are  shown  in 

figures  4,  5,  and  6,  Assuming  a  Gaussian  shape  and  making  use  of  the  fact  that 
the  (accelerated)  molecular  wakes  and  the  turbulent  wakes  are  approximately 

c. 

uncorrelated,  we  obtained  values  of  0^  from  these  curves  in  conjunction 

with  the  mean  temperature  curves.  The  assumptions  involved  in  this  analysis 

i 

*! 

nr*  apparently  not  wholly  satisfied,  since  the  value  of  QC  obtained  depends 

<  _  r 

on  the  particular  points  of  the  correlation  function  used  for  its  computation.  An 

i 

h 

overage  value  of  CP  was  taken  at  each  y  ,  and  these  are  plotted  in 
figure  7,  The  values  of  show  no  consistent  variation  with  y  at  a 

particular  X  ,  and  within  the  scatter  they  are  constant,  indicating  independence 

/ . 

of  the. two  diffusion  phenomena.  The  fact  that  the  ratio 

'i 

unity  for  all  diffusion  times  considered  is  surprising,  and  presumably  a  coincidence 

t 

for  the  particular  experimental  conditions, 

-  ' 

The  mean-wakerat  -  18-7/8  laches  is  shown  in  figure  8  no  typical  of 
these  measurements.  Also  shown  is  a  comparison  of  the  turbulent  wake  as  com- 

'•At* 

putad  from  tho  moan  wake  and  as  measured  directly  with  the  trigger  circuit.  The 

. 

circuit  was  used  here  with  the  signal  clamped  tc  its  lowest  value,  l0e.  the  minimum 

•i 

value  attained  by  tho  signal  over  a  time  interval  of  1/2  second  was  used  as  the 
base  line  voltage.  Near  the  edges  of  the  wake  the  circuit  clamps  on  to  the  ambient 


V 


is  approximately 


temperature  mad  the  two  curves  agree  closely.  When  the  joint  piohnbil  '.ties 
were  measured,  the  circuit  was  adjusted  ao  that  the  output  corresponded  to 
the  correct  turbulent  wake,  aa  explained  in  the  preceding  sec. ion. 

The  joint  probability  functiona  were  obtained  by  measuring  a  set  of  con¬ 
ditional  probability  curves  and  then  using  these  to  construct  a  contour  plot  of 
the  desired  functions.  This  plot  was  normalised  to  unity  at  its  peak  value.  A 
typical  set  of  measured  conditional  probability  curves  is  shown  in  figure  9. 

The  contours  were  constructed  from  interpolation  of  such  data.  The  contours 
were  symmetrical  around  ^  in  all  cases,  as  required  by  isotropy, 

so  feat  only  those  on  one  side  of  y  (  *  were  plotted. 

Near  the  sources,  whsre  the  Lagrangian  velocity  correlation  coefficient  is 

i 

close  to  unity,  the  measured  joint  probability  densities  of  the  displacements 
should  be  identical  with  the  velocity  joint  probability  densities  at  the  two  source 
points  (with  the  appropriate  scaling).  The  measured  joint  probability  densities 
of  the  displacements  at  7/8  inches  (figures  10  and  11)  were  therefore  used  to 
calculate  the  double  correlation  of  the  velocities  fit  the  same  location,  and  these 
values  are  plotted  in  figure  12  together  with  the  directly  measured  velocity 
correlation  curve.  The  computed  points  give  a  smaller  correlation  coefficient 
than  the  direct  meaaurement,  a  circumstance  that  may  bo  attribu.able  to  the 

i  • 

Lagrangian  correlation's  not  being  sufficiently  close  to  one  at  the  X  -location 
used. 

j  Wbsn  molecular  diffusion  has  been  in  action  so  long  that  the  molecular 
wakes  spread  into  each  other  for  an  appreciable  range  of  (y  t  -x)  ,  the 

measuring  technique  fails  since  the  wakes  become  indistinguishable.  In  thi3 

\ 


region  the  joint  probability  (unction  is  not  completely  recoverable .  Because  of 
this  reduced  reliability  the  measurements  made  where  this  effect  was  dominating 
are  not  reported  here.  These  include  measurements  at  V  ■  6-7/8  inches, 
18-7/8  inches,  and  Y  *  1/4  inch. 

The  data  for  Y  ■  3/4  inch  (figures  11, 13, 14)  illustrate  the  propagation 
of  the  Joint  probability  density  for  particle  displacement  in  time.  The  particle 
motions  remained  correlated  for  all  diffusion  times  observed,  as  indicated  by 
the  non-aero  eccentricity  of  the,  elliptical  contours.  At,  X  *  18-7/8  the 
flattening  of  the  contours  near  (  *  Y  is  a  manifestation  of  the  con¬ 

straint  forbidding  particle  paths  to  cross. 

From  the  contour  plots  at  Y  *  3/4",  the  marginal  distributions  of 
were  computed  and  are  shown  in  figure  15.  For  this  particular 
source  spacing  and  at  these  diffusion  times,  this  function  stays  similar  and  is 
approximately  Gaussian  in  shape. 

The  marginal  distribution  of  y  ^  and  ^  computed  from  the 
probability  hills  did  not  correspond  exactly  to  the  directly  measured  distribution 
of  ^  or  .  This  difficulty  is  due  to  the  lack  of  precision  at  low  values 

of  the  probability,  and  possibly  to  a  lack  of  independence  between  the  turbulent 
and  molecular  wakes. 

cv 


The  growth  of  48  *  function  of  X 


is  shown  in  figure  16. 


CONCLUSIONS 


This  study  of  measuring  technique  for  determining  the  joint  probability 
density  of  particle  displacements  in  a  turbulent  flow  field  discloses  several 
effects  that  have  to  be  correctly  accounted  for  in  any  further  work  along  these 
lines.  When  fluid  particles  are  marked  with  a  substance  that  can  diffuse  by 
molecular  motion,  the  resulting  finite  size  of  the  convectcd  region  enters  the 
measurements  in  such  a  fashion  that  the  laminar  and  turbulent  effects  are  in¬ 
separable.  Near  the  sources,  given  y  values  for  the  wake  locations  imply 
a  value  of  the  velocity  derivative.  This  same  velocity  gradient  will  change  the 
widths  of  the  molecular  wakes,  and  this  change  is  reflected  Ln  any  measurements 
weighted  by  the  wake  widths  at  a  given  value  of  the  temperature.  The  large  shift 
of  the  peak  value  of  the  probability  hill  from  the  geometrical  center  (  y  *  if*.  = 
in  the  present  measurements  illustrates  this  effect  most  noticeably.  The  shift  had 
been  entirely  attributed  to  the  known  sktwneas  of  the  velocity  difference  in  isotropic 
turbulence  until  it  was  seen  that  the  marginal  distributiomof  y(  and  y*, 
computed  from  these  data  also  had  the  peaks  shifted.  Assuming  that  the  wakes 
were  acted  on  by  the  same  velocity  gradient  for  all  diffusion  times  corcidnred,  and 
estimating  the  velocity  derivatives  at  each  v/ake,  yields  a  result  io  :  the  peak  shift 
in  rough  agreement  with  the  measured  value  which  indicates  that  the  analysis  has 
a  proper  basis. 

The  results  show  further  more  that  the  average  shape  of  the  molecular 
wake  at  a  given  diffusion  time  la  independent  of  y  ,  if  no  other  condition 
Is  specified,  so  that  the  assumption  of  uncorrclntcd  molecular  and  turbulent  effect3 

Is  usable  in  ths  single  source  case. 


23 

The  observed  changes  in  time  for  the  Y  *  3/4  inch  joint  probability 
densities  show  that  the  Joint  probability  density  function  does  not  stay  similar 
in  time  but  requires  more  than  one  parameter  for  its  description.  This  seems 
reasonable  because  the  constraints  controlling  the  behavior  of  the  function  at 
various  diffusion  times  are  different.  Near  the  source  the  displacement  Joint 
probability  distribution  is  similar  to  the  velocity  joint  probability  distribution, 
while  far  from  the  source  the  fact  that  the  wakes  cannot  cross  becomes  significant. 

i 

.jj  Finally,  it  Is  suggested  that  studies  for  long  diffusion  times  (when  molecular 
broadening  has  "overlapped"  the  two  wakes  in  these  air  experiments)  be  carried 
out  In  water,  which  has  a  considerably  higher  Prandtl  number. 
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